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1 Introduction. 

In this work we consider quadratic Hamiltonians of the form 

H = (M, AM) + (M, Bj) + (7, C7) + (P, M) + (Q, 7). (1) 

where M = {Mi, M2, M^), 7 = (71,72,73) and A,C are constant symmet- 
ric matrices, i? is a general constant matrix and P, Q are constant vectors. 
Without loss of generality we choose 

= diag(ai, 02,03). (2) 
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The equations of motion are given by 



dMi 



{H,M,} 



dli_ 
dt 



dt 



where Poisson brackets are defined by 



{Mi, Mj} = Eijk Mk, {Mi, 7^} = Eijk 7fc, {li, 7j} = (3) 



with £ijk being the totally anti-symmetric tensor. These linear Poisson brack- 
ets are related to the Lie algebra e(3). 

Because they admit the two Casimir functions 



we need only one additional first integral for Liouville integrability ^ . In our 
paper we restrict ourselves to polynomial first integrals. The relevance of the 
combination of quadratic Hamiltonians with Poisson brackets (jH)) arises 
from the Euler-Poinsot model describing motion of a rigid body around a 
fixed point under gravity and from the Kirchhoff model describing the motion 
of a rigid body in ideal fiuid 12 El EI- 

In section |2l we give a list of known integrable quadratic Hamiltonians 
with additional polynomial integrals of motion. Results of our search are 
given in section 01 

It is known from literature that for pairwise non-equal Oi, 02, 03 in (0) no 
extra integrable cases besides those mentioned in section |21 exist. In addition 
all Hamiltonians (0) with linear or quadratic first integral are known too. 
Hamiltonians with cubic first integral are discussed in We therefore 
study first integrals / of fourth degree and consider the case ai = 02 7^ 03. 

The motivation for a systematic investigation of this ansatz arose from 
the recent finding of V. Sokolov of a new integrable Hamiltonian of this form 
with a fourth degree first integral. Our investigation is in so far more gen- 
eral as we consider quadratic Hamiltonians with linear terms. Our complete 
classification of this problem resulted in new integrable Hamiltonians with 
complex coefficients. The simplest example of such a kind is 



which commutes with some fourth degree polynomial under so(3)-brackets. 

We also were inspired by the Goryachev-Chaplygin Hamiltonian [3 E] 
which admits a fourth degree integral exclusively on the Casimir level J2 = 0. 
Unfortunately our complete classification did not result in new integrable 
cases. 



Ji = (7,7), J2 = h,M) 



(4) 



H = Ml + Ml + 2Ml + iMi + M2 



2 



In section HI we follow [T71 HHj and consider the quantum counterpart 

[Mi ,Mj] = £ijk Mk , [Mi , 7j] = £ijk Ik , [li , Ij] =0 (5) 

of the Poisson bracket Here Mj, 7^ are elements of an associative algebra 
with commutator relations (0). The Hamiltonian is a (non-commutative) 
polynomial of second degree and first integrals are polynomials which com- 
mute with the Hamiltonian. Also here we classify quadratic Hamiltonians 
with linear terms having a fourth degree first integral and obeying the restric- 
tion ai = a2 7^ 03. We find four cases, among them the quantum analogue 
of the above mentioned integrable Hamiltonian of V. Sokolov. 

In all three investigations an ansatz for the Hamiltonian and an ansatz for 
a fourth degree polynomial were made, the Poisson bracket or commutator 
were computed and set to zero, resulting in large non-linear over-determined 
algebraic systems. Although they are linear in the coefficients of H and 
linear in the coefficients of the first integral, the algebraic systems are still 
very large (see table 1) and challenging to solve. 

We used the computer algebra program Crack which was originally de- 
signed to solve over-determined PDE-systems. But its interactive capabilities 
allowed to solve large algebraic systems at first interactively and in doing so 
to learn how to take advantage of the bi-linearity. Due to an ongoing effort 
to generalize gathered experience and to incorporate it into the program it 
is now able to solve large bi-linear systems, for example the first and third 
system in table 1 automatically and the second system with only few man- 
ual interactions. An overview of essential features of the program and other 
applications requiring the solution of bi-linear systems is given in [7]. 



type of problem 


e(3) 


e(3),J2 = 


e(3) quant. 


# of unknowns (i^,/,total) 


17,200,217 


17,176,193 


17,200,217 


# of equations 


451 


396 


451 


total # of terms 


5469 


5243 


9681 


average # of terms/equ. 


12.1 


13.2 


21.5 


time to solve 


18h 53min 


15h 


llh 43min 


details in section 








H 



Table 1. An overview of the solved algebraic systems 



Times are measured on a 1.7GHz Pentium 4 running a 120 MByte REDUCE 
session under Linux. 
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2 Known integrable Hamiltonian on e(3) 

In the following we list all known integrable Hamiltonian of type (0) on e(3). 
First let us consider all cases where the matrices A = {aij}, B = {bij} and 
C = {cij} are diagonal, i.e. the Hamiltonian takes the form: 

H = aiMf + asMi + agMf + 2f>iMi7i + 2b2M272 + 263M373+ 

(6) 

Ci7? + C27i + C373 + Pi^i + + P3M3 + gi7i + g272 + ^aTs- 

Kirchhoff 's problem of the motion of a rigid body in ideal fluid. 

In this case the Hamiltonian does not contain linear terms, i.e. we have 
Pi = Qi = 0- For this problem there are three known classical integrable cases 
of Kirchhoff jH!, Clebsch ^ and Steklov-Lyapunov [1011111. 
The Kirchhoff case is defined by the identities 

ai = 02, bi = 62, ci = C2. 

The additional integral / is linear: / = M3. 

For the Clebsch case the coefficients ai are arbitrary and the remaining 
parameters satisfy the following conditions: 

61 = 62 = 

Cl - C2 ^ C3 - Ci ^ C2 - C3 _ ^ 

as a2 ai 

If not all ttj are equal then the Hamiltonian can be represented in the form 

H = aiMl + a2Ml + asM^ + 02037^ + 030172 + 010273. 

The additional quadratic integral admitted by H is 

I = Ml + Ml + Ml + (02 + 03)72 + (03 + 01)7^ + (oi + 02)7|- 

If Oi = 02 = 03 then we get the Neumann Hamiltonian 

H = Ml + Ml + Ml + C172 + C27I + 0373'- 

The additional integral in this case coincides with H . 

For the Steklov-Lyapunov case the coefficients (Xi cirG arbitrary and the 
remaining parameters satisfy the following conditions: 

bi — 62 63 — 61 bo — 63 

+ + = 0, 

O3 O2 Oi 
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Ci = C2 = C3 . 

fli a2 as 

If not all the are equal then the Hamiltonian can be represented in the 
form 

H = aiMf + a^Ml + a^Ml + 2a2a^Mi-ii + 2a^aiM2-i2 + 2aia2M37s 
+ai(a2 - as)^7i + ^2(03 - aifll + ^3(01 - 02)^73- 
The additional integral is quadratic: 

I = Ml + M| + Ml + 2(a2 + as)Mi7i + 2(as + ai)M272 + 2(ai + a2)Ms73 

+ (a2 - 03)^71 + (^3 - ciifll + («i - «2)^73- 

If ai = a2 = cts we have to interchange the Hamiltonian H and the integral 
/ just as in the Clebsch case. 

Recently in the paper ^2] by Sokolov a Hamiltonian with non-diagonal 
matrix B having an integral of fourth degree was presented. One of the 
possible forms of this Hamiltonian is: 

H = Ml + Ml + 2 Ml + 2 (ai7i + 0272) M3 - (a? + al) 73I (7) 

The integral can be represented as a product I = kik2, where ki = M3 and 

k2 = {Ml + Ml + Ml) Ms + 2 (fliMi + a2M2) {M^i + M2I2) 

+2 (ai7i + 0272) Ml + (ai7i + 0272)^ Ms 

-{al + al) (2Mi7i + 2M272 + Ms7s) 73. 

This case appears to be similar in its properties to the Kowalewski case given 
below. 

The problem of motion of a rigid body around a fixed point. 

Hamiltonians describing such situations have the form (jHl), where 6j = q = 
Pi = 0. The following integrable cases are known. 
The Lagrange case: 

H = Ml + Ml + a,Ml + gs73 (8) 

where a^, are arbitrary. The additional integral I is linear: I = M3. 
The Euler case: 

H = aiMl + a2Ml + asM| (9) 
with an additional quadratic integral: 

I = Ml + Ml + Ml. 
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The Kowalewski case: 

H = Ml + Ml + 2Ml + gi7i + ^272 (10) 

with arbitrary parameters gi, q2- The additional integral of fourth degree can 
be written as I = G\ + where 

Gi = Ml - Ml - gi7i + 5272, G2 = 2M1M2 - g27i - gi72- 

We note, that if q2 7^ iqi, then the transformation M T M, T —>■ TT, 
where 

(cos (f) sm(/) \ 
-sin(/) COS0 (11) 
ij 

can be used to make g2 to zero. 

Generalizations. Terms linear in the moments Mj that occur in the 
Hamiltonian could be interpreted as an action of hydrostatic forces (|4j). We 
have no comments on the physical meaning of other mixed Hamiltonians 
(i.e. Hamiltonians having both B ^ or C ^ 0, and Q ^ ). 

For example, an obvious hybrid of Lagrange's and Kirchhoff 's Hamilto- 
nians (with additional hydrostatic member) is 

H = Ml + Ml + siMl + S273M3 + S373 + ^4^3 + S573, (12) 

3 Results 

3.1 The classical case. 

For the Hamiltonian (0) on e(3) we consider the case of A being diagonal: 

A = diag(ai, 02,03), 

where 

ai = a2 ^ as, 7^ 0, ^ = 1, 2, 3. 

Without loss of generality matrix B can be considered as low-triangular, i.e. 

bi2 = bi3 = &23 = 0. 

Note, that the addition of Casimir functions ^ to the Hamiltonian does 
not influence the equation of motion. By subtracting an appropriate linear 
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combination of Ji and J2 we can make bn = cn = and get 
H = M^ + M^ + a3Ml + 

62171^2 + hiliMs + 63272M3 + 62272M2 + 63373^3 + 

C127172 + C137173 + C2272 + C237273 + C3373 + (13) 

PiMi + P2M2 + P3M3 + gi7i + 5272 + g373- 

Lemma. For a Hamiltonian of type (|T3|l . a canonical transformation can 
be used to obtain 621 = 0. 

Thus we will consider only the following Hamiltonian: 

H = + Ml + a^Ml + 63171 M3 + 63272M3 + 62272M2 + 63373M3 + 

C127172 + C137173 + C2272 + C237273 + C3373 + (14) 

Pi Ml + P2M2 + P3M3 + gi7i + §272 + g373- 

Theorem l.The Hamiltonian of the above form commutes with a poly- 
nomial integral of 4*^* degree iff it coincides with one of the following: 

• (O, where St - are arbitrary, or 



H = Ml + M| + M| + 2S173M3 - shI + S27i + S372 + S473+ 
A (2sxMl + S2M1 + S3M2 + S4M3 + si(s27i + S372) j , 

H = Mf + Ml + 2M| + si (z7i + 72) M3+ 

(16) 

S2{-tMi + M2) + S3M3 + S4(^7l + 72) - S1S273, 
H = Mf + Ml + 2M| + 2(si7i + S272) M3 - {sj + sl)-fl+ 

(17) 

S3 (M3 + si7i + S272) + S471 + S572, 
where Si are parameters constrained only by the condition: 

S2S5 + S1S4 = 0. 

Remarks 

1. The case ()15|) was founded by Rubanovskiyi jTHj- 

2. Under the conditions Si = S2 = S3 = the formula (fT7|) describes the 
Kowalewski case. If S3 = S4 = S5 = we have the Sokolov case for which 
the general formula (fT7j) was presented in ^3] and the Lax pair was found 
by Sokolov and A.V. Tsiganov jl5j . 

3. The Hamiltonian (fTHIl is probably new. 
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3.2 Classical case with additional condition. 

We now consider the case of the area integral J2 being equal to zero. In this 
case the following additional integrable Hamiltonians are known. 
The Goryachev-Chaplygin Hamiltonian jSl Ej 

H = Mf + M| + 4M| - si7i (18) 

belongs to the class of Hamiltonians describing the motion of a rigid body 
around a fixed point. The additional integral is of third degree. 

For the Kirchhoff problem is known the Hamiltonian by Chaplygin: 

H = Mf + M| + 2M| + c^j - C7|. (19) 

Here the additional integral has fourth degree. 

Theorem 2. The Hamiltonian of the form (|TH) commutes with a poly- 
nomial integral of fourth degree with additional condition (M, 7) = 0, iff it 
coincides with one of the Hamiltonian from Theorem 1 or with the following: 

• 

H = Ml + Ml + 2M| + 51(7? - 71) + S27172 + S3M3 + S471 + S572, 

(20) 

where Si - are arbitrary parameters. 
The Hamiltonian (j^n|) is the generalization of Chaplygin's case [SI lEl ^] • 

4 Quantum case. 

In the papers by Komarov, Sklyanin ^3 and others quantum gener- 

alizations of classical Hamiltonians on e(3) and so(4) were considered. Sim- 
ilarly, instead of the Poisson bracket (jHl) we investigate the commutator re- 
lations (0) in an associative algebra with generators Mi, M2, M^, 71,72,73- 
We can consider Mi,^i to be operators. Due to commutator relations (0) 
any monomial always could be ordered such, that Mi come before any 7j and 
indices increase (such a presentation is unique). Thus Hamiltonians have the 
same form as in the classical case, but the multiplication is now non- 
commutative. In this case, any element / from the associative algebra that 
satisfies [H, J] = is called an integral. 

Theorem 3. The Hamiltonian of the above form commutes with a poly- 
nomial integral of 4*'* degree iff it coincides with one of the following: 
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H = Mf + M| + siM| + S273M3 + S37I + S4M3 + S573, (21) 
where Si - are arbitrary parameters; 



H = Mf + M| + Ml + 2S173M3 - sIjI + S2I1 + S372 + S473+ 

/ X (22) 

A ( 2siM| + S2M1 + S3M2 + S4M3 + si(s27i + S372) 



(23) 



(24) 



H = Ml + M| + 2M| + si(27i + 72)M3+ 

S2(-iMi + M2) + S3M3 + 54(^71 + 72) - S1S273, 

H = Mf + M| + 2M| + 2(si7i + S272) M3 - (s? + ^2)72+ 

S3 (M3 + si7i + S272) + S471 + S572, 
where Si are parameters constrained only by the condition: 

+ S1S4 = 0; 

Notice that this hst looks exactly as the list of subsection 13.11 although a 
priori this coincidence is not predictable. The most interesting here is Hamil- 
tonian (j2H), which is a quantum generalization of the Sokolov Hamiltonian. 
If Si = S2 = S3 = we have the quantum Kowalewski Hamiltonian [T71 l2Uj . 



Acknowledgements 

The authors are grateful to V.V. Sokolov for useful discussions. The research 
was partially supported by RFBR grant 02-01-06670. 



References 

[1] Arnold V.I. Mathematical methods of classical mechanics, Moscow: Sci- 
ence, 1974, 431 p.. 

[2] Kozlov V.V. Symmetries, topology and resonances in Hamiltonian me- 
chanics, Publ. Udmurdski state university: Izhevsk, 1995. 



9 



[3] Borisov A.V., Mamaev I.S. Rigid body dynamics, Publ. Udmurdski state 
university: Izhevsk, 2001. 

[4] Borisov A.V., Mamaev I.S. Poisson's structures and algebra Lie in 
Hamiltonian mechanics, Udmurdski state university: Izhevsk , 2001, 460 
P- 

[5] Goryachev D.N. New cases of motion of a rigid body around a fixed 
point , 1915, book 3, p. 1-11. 

[6] Chaplygin S.A. A new particular solution of the problem of motion of 
a rigid body in fluid, Sobranie sochineny v. 1,M.-L.: GITTL, 1948, 
p.337-346. 

[7] Wolf, T. Applications of Crack in the Classification of Inte- 
grable Systems, preprint, submitted to CRM Proceedings, 15 pages, 
nlin.SI/0301032 (2002). 

[8] KirchhofT G.R. Vorlesungen iiber mathematische Physik. Mechanik, 
Leipzig, 1874. 

[9] Clebsch A. Uber die Bewegung eines Korpers in einer Fliissigkeit, Math. 
Annalen, Bd.3, 1871, s.238-262. 

[10] Steklov V.A. On the motion of a rigid body in a fluid, Kharkov, 1893, 
234 p.. 

[11] Lyapunov A.M. A new case of integr ability for the equation of motion 
of a heavy rigid body in fluid, Sobranie sochineny v. 1, M., 1954, p. 
320-324. 

[12] Sokolov V.V. A new integrable case for the Kirchhoff equation, Teoret. 
and Mat. Fiz., 2001, 129, no.l, 31-37, translation in Theoret. and Math. 
Phys., 2001, 129, no.l, 1335-1340. 

[13] Rubanovskiyi V.N. New cases of integrability for the equation of motion 
of a heavy rigid body in fluid, Vestnik MGU, ser.math. mach., 1968, 2, p. 
99-106. 

[14] Sokolov V.V. A generalized Kowalcwski Hamiltonian and new integrable 
cases on e(3) and so(4), in the book "Kowalevski property", edt. V.B. 
Kuznetsov, to appear in CRM Proceedings and Lecture Notes, AMS, 
2002. Preprint nlin.SI/0110022, 2001. 



10 



[15] Sokolov V.V., Tsiganov A.V. On the Lax pairs for the generahzed 
Kowalewski and Goryachev-Chaplygin tops, Teor.Math.Phys., v. 131(1), 
p. 543-549, 2002. 

[16] Yehia H.M. New integrable problems in the dynamics of rigid bodies 
with the Kovalevskaya configuration. I- The case of axisymmetric forces, 
Mech. Res. Com., 1996, v.23, 5, p.423-427. 

[17] KomarovI.V. Kovalevskaya basis for atom of hydrogen, Teor.Math.Phys, 
1981, V.47, p.67-72. 

[18] Komarov l.V. Goryachev-Chaplygin top in quantum mechanics, Teor. 
Math. Phys, 1982, v.50, p.402-409. 

[19] Sklyanin E.K. Goryachev-Chaplygin top and quantum method of re- 
verse scattering task, in the book "Differential geometry, Lie groups and 
mechanics" L.: Science, 1984, p. 236-257 

[20] O.Lapportc Note on Kowalewski's top in quantum mechanics. Phys. 
Rev., 1933, v.43, p.548-51 



11 



